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A.  PURPOEC  OF  PROGRAM 
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To  compute  the  specific  energy  and  specific  force  in  a  circular  open  channel. 
Reference:  Chow,  Ven  Te,  Open  Channel  Hydraulics,  McGraw-Hill,  1959 »  PP  ^1,5^* 


«.  PROGRAM  SPECIFIC  AT  IONS 
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c.  methods 

The  program  is  written  in  G635  time-share  series,  FORTRAN  IV,  and  is  part  of 
a  Conversationally  Oriented  Real-Time  Program-Generating  System  (CORPS).  It 
consists  of  a  main  program  and  a  subroutine.  All  I/O  requirements  are  handled 
in  the  main  program  with  the  subroutine  handling  the  computation. 
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equipment  details 

The  program  was  developed  and  is  operational  on  WES  G635>  Vicksburg,  MS.  It 
is  also  operational  on  the  HIS  66/80,  Macon,  GA,  and  Boeing  CDC ,  Seattle,  VA. 


"irndwt-60'f»uT - - - - 

The  required  inputs  which  are  entered  at  execute  time  are:  discharge  cfs, 
diameter  ft,  number  of  flow  depths,  and  flow  depths,  ft.  The  output  in¬ 
cludes  the  given  data,  plus  the  specific  energy  and  specific  force. 
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Complete  documentation  of  this  program  is  available  from  the  Engineer 


PROGRAM  SPEC! F 1 CAT I  ON 


LariKurij'e :  ANSI  FORTRAN  (FORTRAN  IV) 

Solution  Requirements :  The  run  command 

RUN  WEGL1 B/CORPS/H6125 ,  R 


plus  the  input  variables  defined  in  (E). 

Method  of  Analysis:  Solves  algebraic  equations  for  the  specific- 
energy  and  specific  force. 

Core  Requirements  G635:  10  K  words 

External  Storage:  None 

Restrictions :  Velocity  distribution  is  uniform  across  the  flow 
section  and  the  invert  slope  is  small  (i.e.,  <  10  degrees). 
General  Equations: 


Specific  Energy  E  (ft  of  water) 


E 


y 


O  -  0 

Specific  Force  F  (ft  of  water)  F  =  Ay  +  — 

where:  Q  is  the  discharge  (cfs),  g  is  acceleration  of  gravity 
(32.2  ft/sec  ),  y  is  flow  depth  (ft),  A  is  cross-section  area 
(ft^)  and  y  is  the  distance  (ft)  of  the  centroid  of  the  area  A 
below  the  surface  of  flow. 

Range  of  Quantities:  Unlimited  for  practical  application. 

Accuracy:  Governed  by  accuracy  of  input  data;  the  specific  energy 
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and  specific  force  are  printed  to  +0.01  ft  and  ft  . 
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HEF:  EH  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  Computer 

Program  Library  Standards  and  Documentation,  Appendix  B  V 

PART  I:  ENGINEERING  DESCRIPTION 

1.  PROGRAM  NUMBER:  722-F3-RO-6AK 

2.  TITLE :  H6125  -  Specific  Energy  and  Specific  Force  in  a  Circular 

Open  Channel. 

3.  REVISION  LOG:  N/A 

^ .  PURPOSE  OF  PROGRAM:  To  compute  the  specific  energy  and  specific 
force  in  a  circular  open  channel. 

Reference : 

Chow,  Ven  Te ,  Open  Channel  Hydraulics,  McGraw-Hill,  1959 9  PP  ^1,5^* 
5.  STEP  SOLUTION: 

a.  The  inputs,  discharge  (Q)  cfs,  diameter  (D)  ft,  and  flow  depth 
(y)  ft,  are  entered. 

b.  The  computational  steps  are: 

(1)  Area  (A)  ft2 


(2)  Moment  (Ay)  ft^  of  area  with  respect  to  the  surface 


(3)  Specific  Energy  (E)  ft  of  water 

E  =  y  +  £  ■■■ 

2gA2 

2 

where  g  is  acceleration  of  gravity  (32.2  ft /sec  ) 
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(4)  Specific  Force  (F)  ft  of  water 

- 

F  =  Ay*S_ 

c.  The  given  inputs,  plus  the  specific  energy  and  specific  force 
are  printed. 

6.  ACCURACY :  Governed  by  input  data;  the  specific  energy  and  specific 
force  are  printed  to  the  nearest  hundredth  of  a  foot  and  hundredth  cubic 
foot. 

7.  REMARKS:  Velocity  distribution  is  uniform  across  the  flow  section 
and  the  invert  slope  is  small  (i.e.,  <  10  degrees). 
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3-  EQUIPMENT  AND  OPERATING  SYSTEM:  The  program  was  developed  on  a 
G635  time-share  system  in  which  input/output  equipment  consisted  of  a 
Model  33  remote  teletype.  It  is  now  operational  on  the  WES  G635» 
Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 

U.  INPUT  REQUIREMENTS:  The  required  inputs  are  entered  via  the  user’s 
time-share  terminal  device  in  free  field  format.  All  input  cues  and 
reads  are  performed  in  the  main  program.  The  subroutine  handles  the 
computation.  Since  computations  a;  done  in  subroutine  H6125,  the 
necessary  inputs  to  the  subroutine  are  passed  via  the  CALL  statement. 

The  calling  sequence  is: 

CALL  H6125  (arg^, . .  .  ,arg<- ) 

where : 

arg^  -  discharge,  cfs 

arg^  -  diameter,  ft 

arg^  -  flow  depth,  ft 

arg^  -  specific  energy,  ft  of  water 
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arg^  -  specific  force,  ft  of  water 
All  arguments  are  floating  point.  Arguments  1-3  are  the  required  in¬ 
puts  and  arguments  U  and  5  are  the  outputs. 

5.  SECONDARY  STORAGE  INPUT  FORMAT:  None 

6.  INPUT  DATA  DESCRIPTION:  The  following  names  are  used  for  the  input 
variables  in  program  H6l25* 

DISCH  -  discharge,  cfs  (arg^ 


k 


DIAMTR  -  diameter,  ft 


LQZ 


-  equal  1,  execute  all  input  cues  and  reads;  equal  2,  call 
WESLIB  routine  RERUN  and  enter  only  desired  inputs 

LQX  -  equal  1,  print  instructions  from  RERUN;  equal  3,  no  print 

JKL  -  direct  return  from  RERUN  to  desired  input  read 

KKK  -  total  number  of  inputs;  passed  to  RERUN 

WK1  -  working  storage,  one-half  top  width,  ft 

WK2  -  working  storage,  ft 

WK3  -  working  storage 

ZZZZZ  -  character,  equal  RE,  rerun  program;  equal  ST,  stop 

10.  EXAMPLE  CASE:  Compute  the  specific  energy  and  specific  force  for 
6  given  flow  depths . 
a.  Input  data: 

Discharge  (DISCH)  =  100.00  cfs 
Diameter  (DIAMTR)  =  5 <00  ft 
Number  of  flow  depths  (N)  =6 

Flow  depths  (Y(l),  for  I  =  1,  N)  =  1.00,  1.1*3,  2.50,  4.00,  4.50, 
and  5 <00  ft 


b .  Output : 


H6l2 
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INPUT  H6I25-SPECIFIC  ENERGY  AND  SPECIFIC  FORCE  IN  A  CIRCULAR  OPEN 
CHANNEL 

AA-ENTER  DI SCHARGE , CFS . 

=  100 

AB-ENTER  DIAMETER, FT. 

=  5 

AC-ENTER  THE  NUMBER  OF  DEPTH ( S )  FOR  WHICH  SPECIFIC  ENERGY  AND  FORCE  ARE 
TO  BE  CALCULATED.  MUST  NOT  EXCEED  25  DEPTHS. 

=  6 

AD-ENTER  THE  6  DEPTH ( S )  SEPARATED  BY  COMMAS. 

=1,1.43,2.5,4,4.5,5 


OUTPUT  H6125-SPECIFIC  ENERGY  AND  SPECIFIC  FORCE  IN  A  CIRCULAR  OPEN 
CHANNEL 


DISCHARGE  = 

100  . 

.00  CFS 

DIAMETER  = 

5. 

.  00  FT 

FLOW 

SPECIFIC 

SPECIFU 

DEPTH 

ENERGY 

FORCE 

:ft) 

(FT) 

( FT  *  #3 ) 

1  .00 

20.87 

112.23 

1  .43 

8.66 

69.72 

2.50 

4.11 

42.05 

4.00 

4.55 

49.03 

4.50 

4.95 

56 . 16 

5.00 

5.40 

64.90 

ENTER  RERUN 

OR  STOP 

=  STOP 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  and  Computer 
Program  Library  Standards  and  Documentation,  Appendix  C 

PART  Ills  FILE  DOCUMENTATION 

1.  REVISION  LOG:  N/A 

2.  TITLE:  H6l25  -  Specific  Energy  and  Specific  Force  in  a  Circular 
Open  Channel 

3.  SOURCE  PROGRAM  LISTINGS:  See  pages  9-11 

h.  NUMERICAL  AND  LOGICAL  ANALYSIS:  No  special  numerical  techniques 
used. 

5 .  SUBROUTINES  NOT  DOCUMENTED  IN  ABSTRACT:  None 

6.  MISCELLANEOUS :  The  program  is  part  of  the  CORPS  computer  system. 
CORPS  is  an  acronym  standing  for  Conversationally  Oriented  Real-Time 
Program-Generating  System.  The  program  is  now  operational  on  the  WES 
G635,  Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 
The  source  listing  on  page  9  contains  the  first  line  run  command  and 
brief  for  H6125-  This  first  line  run  command  runs  the  binary  H6125B  of 
the  source  listing  on  pages  10-11  (Fortran  source  of  H6125)  and  attaches 
the  WESLIB  routine  HACCT  and  RERUN. 
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0001* #RUN  WESIIB/C0RP5/'H6125B,R;WE5LIB/RERUH,R;WE5LIB/HACCT,R 

0300  6 1  THIS  PROGRAM  COMPUTES  SPECIFIC  ENERGY  AND  SPECIFIC  FORCE  IN  A 

0805  23  CIRCULAR  OPEN  CHANNEL. 

0310  36VELOCITY  DISTRUBUTION  COEF.  =  UNITY. 

0815  56INPUTS  REQUIRED  ARE  DISCHARGE-CFS ,  DIAMETER-FT,  AND  FLOW 
0820  09DEPTH-FT. 

0825  5 10UTPUT  INCLUDES  THE  GIVEN  DATA  AND  SPECIFIC  ENERGY 
0830  19AND  SPECIFIC  FORCE. 

0999X06FINISH 
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0  0  0  0  1  # tRUN  #= ;/C0RPS/H6125B(N0G0) 

1  0  0  00  CHARAC T ER#5  HFILE 
10010  DIMENSION  Y(25> 

10020  HFILE=5MH6125 

10030  LQZ-1 ; IQX-l 

10040  15000  CALL  HACCT (HFILE) 

10C50  PRINT  15111 

1  006  0  151  1  1  FORMAT  (/’’INPUT  H6 1 25-SPEC  I FIC  ENERGY  AND  SPECIFIC  FORCE  IN 
1  0  07  0*  A  CIRCULAR  OPEN"/ 1  4X  ,  '’CHANNEL  "//  ) 

10080  GO  TOC150O3, 15016), LQZ 
10090  15003  PRINT  15004 

10100  15004  FORMAT ("AA -ENTER  DI SCHARGE , CFS . ” ) 

10110  15005  READ , DI SCH 

10120  GO  TOC  15007 , 150  16  ) ,LQZ 

10130  15007  PRINT, "AB-ENTER  DIAMETER , FT . " 

10  14  0  1  5  0  OB  READ  »  DI AMT  R 
1C150  GO  TO( 15010, 15016) , LQZ 
10160  15010  PRINT  15121 

10165  15121  FORMAT ( "AC- ENTER  THE  NUMBER  OF  DEPTH ( S )  FOR  WHICH  SPECIFIC  E 
1 0 167 INERGY  AND  FORCE  ARE”/” TO  BE  CALCULATED .  MUST  NOT  EXCEED  25  DEPTHS. 
10169*") 

10170  15011  READ.H 

10180  GO  TOC  15013, 15016), LQZ 

10190  15013  PRINT  15113, N 

10195  15113  FORMAT (” AD- ENTER  THE  ”,I2,”  DEPTH  C  5 )  SEPARATED  BY  COMMAS.”) 

1  0200  1  50  14  READ  *  C  Y ( I ) , I  =  1 , N ) 

10210  GO  TOC15019, 15016), LQZ 

10220  15016  KKK  =  4 

10230  CALL  RERUN ( KKK , LQX, JKL ) 

10240  GO  TOU  5005, 1  5003, 1501  1, 1  5014, 1501  9), JKL 
10250  15019  PRINT  1 50 2 0 , DISCH , DI AMIR 

1026  0  15020  FORMAT (/"OUTPUT  H6 1 2 5 - S PEC  I F IC  ENERGY  AND  SPECIFIC  FORCE  IN 
10270*  A  CIRCULAR  OP EN"/ 1 5X , "CM ANN EL "//"DI SCH ARGE  =  ”,F12.2,”  CFS"/” 

1 0275*DI AMETER  = 

10  280*  " , F 1 3 . 2 , "  F T"//4X , " FLOW” ,4X, "SPECIFIC", 3X,"5PECIFIC"/4X, "DEPTH", 4 
1  0 290*X, "ENERGY”, 5X , ”FORC E”/4X , ” (FT)”, 6X.”C FT ) ” , 5X , ” C FT # # 3 ) ” ) 

10310  DO  15027  1=1, N 
10320  DEPT  f I  =  Y C I ) 

10330  CALL  H6125CDISCH, DIAMTR, DEPTH, ENERGY, FORCE) 

10340  PRINT  15028, DEPTH, ENERGY, FORCE 

10350  15027  CONTINUE 

10360  15028  FORMAT (F8.2,2Fll .2) 

10370  PRINT,”  ” 

10380  LQZ  =  2 

10390  CHARACTER  ZZZZ ZZ#2 

10400  16000  PRINT,  "ENTER  RERUN  OR  STOP” 

10410  READ  16001,  ZZZZZZ 
10420  16001  FORMAT ( A2 ) 

10430  I FCZZZZZZ . EO . 2HRE)  GO  TO  15000 
10440  IFCZZZZZZ.EQ. 2H5T)  GO  TO  20000 
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10450  PRINT, "ERROR  **#  RETYPE” 

10460  GO  TO  16000 
10470  20000  STOPiEND 

30  0  0  0  SUBROUTINE  H6 125( DISCH , DIAMTR , DEPTH , ENERGY, FORCE  ) 

30  0  1  0  G-32.2;WKl=:,QRT(DEPTH*DIAmR-DEPTH*#2) 

3002  0  WK7  =  2 . *DCP  TU-DI AMI R ; UK3 - ARCOS< -WK2/DI AMTR ) 

30  0  30  A  =!*JK  1  *k'K2/2  .  +WK 3* DIAMTR* *2/4  . 

30  04  0  AY =WK1 «<  DEPTH *#2/3. -DE P  HI* DI AMTR/3 .  +DI AMTR#*2/4  .  ) *DI AMTR**2*UJK2*U 
30  050*K3/fi  . 

30060  ENERGY-DEPTH  +DI 5CH#*2/( 2 . #G#A##2 ) 

30070  FURCE=AY*DI5CH##2/(G*A) 

300SQ  RETURN 
30090  END 
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